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Abstract 

We consider a two-chain, spin-^ antiferromagnetic Heisenberg spin ladder in an external 
magnetic field H. The spin ladder is known to undergo second order quantum phase transi- 
tions (QPTs) at two critical values, Hci and Hc2, of the magnetic field. There are now known 
examples of strongly coupled (rung exchange interaction much stronger than nearest-neighbour 
intrachain exchange interaction) organic ladder compounds in which QPTs have been exper- 
imentally observed. In this paper, we investigate whether well-known bipartite entanglement 
meaures like one-site von Neumann entropy, two-site von Neumann entropy and concurrence 
develop special features close to the quantum critical points. As suggested by an earlier theo- 
rem, the first derivatives of the measures with respect to magnetic field are expected to diverge 
as — > Hci and H Hc2- Based on numerical diagonalization data and a mapping of the 
strongly-coupled ladder Hamiltonian onto the XXZ chain Hamiltonian, for which several an- 
alytical results are known, we find that the derivatives of the entanglement measures diverge 
as -ff ^ Hc2 but remain finite as H ^ Hd. The reason for this discrepancy is analysed. We 
further calculate two recently proposed quantum information theoretic measures, the reduced 
fidelity and reduced fidelity susceptibility, and show that these measures provide appropriate 
signatures of the QPTs occuring at the critical points H = Hci and H = Hc2- 

PACS number(s): 03.67.Mn 



I. INTRODUCTION 

Antiferromagnetic ( AFM) Heisenberg ladders are examples of interacting many body systems which 
exhibit a range of novel phenomena [11 [2] . An n-chain spin ladder consists of n chains coupled by 
rungs, the simplest example being a two-chain ladder with n — 2. The study of ladders as prototypi- 
cal many body systems became important after the discovery of high temperature superconductivity 
in the strongly correlated cuprate materials. The dominant electronic and magnetic properties of 
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the cuprates are associated with the Cu02 planes which have the structure of a square lattice 
[3]. The level of rigour that can be achieved in the treatment of strong correlation is less in two 
dimensions (2d) than in Id. Ladder models, with structure interpolating between Id and 2d, serve 
as ideal candidates to address issues related to strong correlation and also to investigate how elec- 
tronic and magnetic properties change as one progresses from the chain to the plane. In undoped 
ladder models, each site of the ladder is occupied by a spin (usually of magnitude 5) and the spins 
interact via the AFM Heisenberg exchange interaction. In doped ladder models, some of the spins 
are replaced by positively charged holes which are mobile. The Hamiltonian describing the doped 
systems are the t-J and Hubbard ladder models [DEI [3]. With the discovery of a large number of 
materials having a ladder-like structure, the study of ladders has acquired considerable importance. 
The materials exhibit a range of phenomena including superconductivity in hole-doped systems, 
the 'odd-even' effect in which the excitation spectrum of an n-chain ladder is gapped (gapless) if n 
is even (odd) and quantum phase transitions (QPTs) tuned by an external magnetic field [Tl[2l[3]. 
Many of the experimental observations were motivated by theoretical predictions, superconductivity 
being a prime example [11 IH . 

In this paper, we focus on QPTs in a two-leg AFM Heisenberg ladder (Fig. 1) in an external 
magnetic field. The Hamiltonian describing the model is given by 

L L 

n = 5][j||(Si^,.Si,,+i + S2,,.S2,,+i) + JiSi^,.S2,,] -HY^isi^ + si^) (1) 

where the indices 1 and 2 distinguish the lower and upper legs of the ladder and j labels the rungs. 
The spins have magnitude \ (|Sj| — \) and interact via the AFM Heisenberg exchange interaction. 
The intrachain and rung exchange couplings are of strengths J\\ and Jj. respectively. The total 
number of rungs is L and periodic boundary conditions are assumed. The factor g/is (<7 is the 
Lande splitting factor and /is the Bohr magneton) is absorbed in H . If = 0, the ladder decouples 
into two non-interacting spin-^ Heisenberg chains with no gap to spin excitations. For any arbitrary 
J_L^ 0, the excitation spectrum acquires a gap (spin gap). In the strong coupling limit, Jj. >> J||, 
a simple physical picture of the ground state and the origin of the spin gap can be given. The spins 
along the rungs predominantly form singlets in the ground state. A spin excitation is created by 
replacing a singlet by a triplet which propagates along the ladder due to the intrachain exchange 
interaction. In first order perturbation theory, the spin gap A is given by A w Jj^ — J|| separating 
the lowest excited state from the dimerized ground state. 

There are now several known strong coupling ladder compounds [Tj. Of these, the organic 
ladder compounds Cu2{CzHi2N2)2Ch [3, {C5Hi2N)2CuBri [9] and {bIAP)2CuBri.2H20 [lO] 
are of special interest because of the experimental observation of the QPTs in these systems by the 
tuning of an external magnetic field. A QPT occurs at T = and brings about a qualitative change 
in the ground state of an interacting many body system at a specific value gc of the tuning parameter 
9 [H]- QPTs are driven by quantum fiuctuations and in the case of second order transitions, the 
quantum critical point is associated with scale invariance and a diverging correlation length. The 
ground state energy becomes non-analytic at the critical value Qc of the tuning parameter. If one 
of the phases is gapped, the gap goes to zero in a power-law fashion as, g ^ gc- In the case of 
a spin ladder, the external magnetic field H plays the role of the tuning parameter g. There are 
two critical points, Hci and Hc2 [IllHll9l[I^. AtT = and for < < Hci, the ladder is in 
the spin gap phase. In the presence of the magnetic field, there is a Zeeman splitting of the triplet 
{S = 1) excitation spectrum with the S"^ = 1 component having the lowest energy. The spin gap 
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is now A — H. At H = Hd = A, the gap closes and a QPT occurs to the Luttinger liquid (LL) 
phase characterized by a gapless excitation spectrum. At the upper critical field H — Hc2, there 
is another QPT to the fully polarized ferromagnetic (FM) state. The magnetization data exhibit 
universal scaHng behaviour in the vicinity of Hd and Hc2, consistent with theoretical predictions 
[HIHllSlIini- In the gapless regime Hd < H < Hc2, the ladder Hamiltonian can be mapped onto an 
XXZ chain Hamiltonian the thermodynamic propeties of which can be calculated exactly using the 
Bethe Ansatz (BA) [U [U [9l [10] . The theoretically computed magnetization versus magnetic field 
curve is in excellent agreement with the experimental data. QPTs can be observed in the organic 
ladder compounds as the magnitudes of the critical fields are experimentally accessible. 

In recent years, QPTs have been extensively studied in spin systems using well-known quantum 
information theoretic measures. A number of entanglement measures have been identified which 
develop special features close to the quantum critical point |13[ [T4j [TSj [T6j \T7\ [TSt I19j . It has 
been shown [16] that, in general, a first order QPT linked to a discontinuity in the first derivative 
of the ground state energy, is signalled by a discontinuity in a bipartite entanglement measure 
whereas a discontinuity or a divergence in the first derivative of the entanglement measure marks a 
second-order phase transition characterized by a discontinuity/divergence in the second derivative 
of the ground state energy. Another measure which provides a signature of QPTs is that of ground 
state fidelity |20[ I21j. The utility of the measure and a related measure, fidelity susceptibility, has 
been explored in a number of studies [221 EH [25l ESI EZl EH] ■ Fidelity, a concept borrowed from 
quantum information theory, is defined as the overlap modulus between ground states corresponding 
to slightly different Hamiltonian parameters. The fidelity typically drops in an abrupt manner at a 
critical point indicating a dramatic change in the nature of the ground state wave function. This is 
accompanied by a divergence of the fidelity susceptibility. In these approaches, the fidelity measure 
involves global ground states. Recently, the concept reduced fidelity (RF) (also called partial state 
fidelity) has been developed, which relates to the fidelity of a subsystem [25l ESI [30l [3ll [32l |33] , along 
with the associated notion of reduced fidelity susceptibility (RFS). Using the RF and RFS measures, 
QPTs have been studied in spin models like the Lipkin-Meshkov-Glick model |30[ [31] , the transverse 
field Ising model in Id [32] and the spin-^ dimerized Heisenberg chains [33]- In this paper, we use 
some well-known bipartite entanglement measures, which include one-site entanglement, two-site 
entanglement and concurrence, for the study of QPTs in the S = ^ two-leg AFM Heisenberg ladder 
(Fig. 1) described by the Hamiltonian given in Eq. (1). We show that the entanglement measures 
develop characteristic features close to the quantum critical point H = Hc2 but not at the critical 
point H — Hd- We next show that the measures based on the RF and RFS signal the occurence 
of QPTs at both the critical points H — Hd and Hc2- 

II. ENTANGLEMENT AND FIDELITY MEASURES PROB- 
ING QPTS 

We first define the various entanglement and fidelity measures which provide the basis of our 
calculations. The single-site von Neumann entropy, a measure of the entanglement of a single spin 
with the rest of the system, is given by 

Sii) = -Trp{i)log2p{i) (2) 

where p{i) is the single-site reduced density matrix [lH[l9]. The two-site entanglement S{i,j) is a 
measure of the entanglement between two separate spins, at sites i and j, and the rest of the spins 
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|181 119j. Let p{i,j) be the reduced density matrix for the two spins, obtained from the full density 
matrix by tracing out the spins other than the ones at sites i and j. The two-site entanglement is 
given by the von Neumann entropy 

S{i,j) = ~Trp{i,j)log2p{iJ) (3) 

In a translationally invariant system, S depends only on the distance n — — A knowledge of the 
two-site reduced density matrix enables one to calculate concurrence, a measure of entanglement 
between two spins at sites i and j [34l|35]. Let p{i,j) be defined as a matrix in the standard basis 
{ITT) ) ITi) ; UT) lii)}- One can define the spin-reversed density matrix asp — [ay®ay)p* [ay®ay), 
where ay is the Pauli matrix. The concurrence C is given by C = max{\i — A2 — A3 — A4, 0} where 
Ai's are square roots of the eigenvalues of the matrix pp in descending order. C = implies an 
unentangled state whereas C — 1 corresponds to maximum entanglement. 

The fidelity F is given by the modulus of the overlap of normaHzed ground state wave functions 
|V'o{A)) and |V'o(A -I- 5X)) for closely spaced Hamiltonian parameter values A and A-|-i5A [20 1 [211 [22] , 

F(A, A + a) = |(V'o(A)|^o(A + (5A))| (4) 

Eq. (4) gives a definition of the global fidelity. The reduced fidelity (RF) [Ml [30l EH [33] refers 
to a subsystem and is defined to be the overlap between the reduced density matrices p = p{h) and 
p = p(h + 5) of the ground states \<j>Q{h)) and \<j>o{h + S)), h and h + 5 being two closely spaced 
Hamiltonian parameter values. The RF is 

FR{h, h + 5)=TryJpippi (5) 

We now compute the different entanglement and fidelity measures for the ladder Hamiltonian 
given in Eq. (1). The external magnetic field H serves as the Hamiltonian parameter. One 
notes that the z-component, S^^i = + '^i.j)^ of the total spin is a conserved quantity. 

Using this fact, the Hamiltonian is diagonalized for different values of L with the help of the 
numerical diagonalization package TITPACK [36]. We take J± ~ 13 K and J|| = 1.15 K which 
are the approximate values of the rung and intrachain exchange couplings in the AFM compound 
[hi AP)2CuBri.2H20 [10]. We determine the ground state as well as the three lowest excited 
state energies for different values of H with L ranging from L = 2 to L = 16. Using the data, we 
examine the variation of the fidelity F{H, H + 5) (Eq. (4)) with increasing magnetic field strength 
H and 6 — .001. We observe sharp drops in F{H, H + 5) at H^^ = Al (inset of Fig. 2), where 
Al is the spin gap, i.e., the difference in the energies of the first excited and ground states. A 
polynomial fitting of the Al versus ^ data points yields « 11.8416 -t- .9739(i) + .6621(-i)2. 
the thermodynamic limit L 00, the critical field is thus Hd — Aqo ~ 11.8416. In the case of the 
strong coupling ladder (Jj_ >> J||), the critical field Hd — J± — J\\ to the first order in perturbation 
theory [37J. The fully polarized FM ground state {H > Hc2) becomes unstable when the lowest 
energy of the spin waves falls below the energy of the polarized state. The magnitude of Hc2 can be 
calculated exactly as Hc2 = J± -|- 2 J| | . The estimates of Hd and Hc2 are in close agreement with 
the experimental results [2 [3 [37]. The numerical diagonalization data reproduces the exact value 
of Hc2 in the thermodynamic limit. We further obtain the variation of the magnetization m{H) and 
its first derivative jjj with H in the thermodynamic limit adopting the extrapolation procedures 
outlined in [38j. The magnetization m{H) is the average magnetization per site and because of 
translational invariance m{H) = {Sf). At T = 0, the expectation value is calculated in the ground 
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state. The inset of Fig. (3) shows that the derivative ^ tends to diverge as H ^ Hd and Hc2- 
This is consistent with the existence of a square root singularity in m{H) in the vicinity of the 
quantum critical points Hd and Hc2 O [37] . Since second order QPTs occur at the critical fields 
Hci and Hc2, the first derivatives of the entanglement measures, S{i), S{i,j) and C, with respect 
to the tuning parameter H may exhibit a discontinuity or a divergence as the critical points are 
approached [16]. We compute the various first derivatives to ascertain that this feature of critical 
point transitions holds true in the case of the spin ladder. 

The single-site reduced density matrix p{i) can be written in terms of the spin expectation value 

{sn as m 

pi^)-( ^""i^^^^ ) (6) 





in the |t) , U) basis. From Eq. (2), 

S{t) = - ^ A, log2 X^ (7) 

i 

where the Ai 's are the two diagonal elements of p{i). Fig. 3 shows the variation of with H. It 
is observed that unlike ^ , ^^j^ tends to diverge only near Hc2, while it approaches a finite value 
close to Hci- The values of Hd and Hc2 are Hd = 11.8416 K and Hc2 = 15.3 K as obtained from 
numerical diagonaHzation data. The strongly coupled ladder model in high magnetic field can be 
mapped onto a Id XX Z AFM Heisenberg chain with an effective Hamiltonian [37l [39] 

L L 

Heff = J|| Et^i + ^'^+1 + 2^^^Ui - ^1 (8) 

where H = H — J± — ^ is an effective magnetic field and 5" 's {a = x, y, z) are pseudo spin- 
i operators which can be expressed in terms of the original spin operators. There are several 
exact results known for the XX Z spin-^ chain in a magnetic field [40^ 141). In particular, the zero 
temperature magnetization m{H) close to the quantum critical points is given by the expressions 
(we use the symbol H instead of H) 

/2 

m{H) ^ ^^{H - Hd)lJ\\, H > Hd (9) 

V2 



m{H)^l-^J{H,2-H)/J\\, H <H,2 (10) 

Similar expressions are obtained in the case of an integrable spin ladder model with the help of 
the thermodynamic BA f^. Using the analytic expressions of m{H) in Eqs. (9) and (10), the 
first derivative of single-site von Neumann entropy with respect to magnetic field H, ^§jf-, can be 
calculated analytically from Eqs. (2) and (6). Again, the derivative diverges near Hc2 (Fig. 4) but 
not as the quantum critical point Hd is approached, consistent with numerical results. The values 
of Hd and Hc2 are Hd = J± — J\\ = 11.85 K and Hc2 = Jj_ -I- 2 J|| = 15.3 K. The estimate of Hd 
is from first-order perturbation theory. 

The correlation functions of the S=^ XXZ chain in a magnetic field are known [42] in the 
gapless phase Hd < H < Hc2- In terms of the original spin operators, these are given by 



5 



1 / 1 \ 

{S!{r)S!{Q)) = — + + cos{27rmr) ( - j (11) 



2K+1 



(S'+(r)5f (0)) cos[tt{1 - 2m)r] (^^^ + cosinr) (J-^ (12) 

where K is the LL exponent. For simpHcity, we have dropped some prefactors (constants) in the 
terms appearing in Eqs (11) and (12). The expressions for the correlation functions are utilized to 
study the variation of the two-site entanglement S{i,j) and concurrence Ci,i+i with respect to the 
magnetic field. These quantities can be computed from the two-site reduced density matrix p{i,j) 
which, in terms of the spin expectation values and correlation functions, is given by [43] 



p{hi) = 



( \ + {SD + (5f 5|) \ 






\ - {SlSf) (SfSj) + (SfS]) 

V i - (5f) + (SfS^) J 

S{i,j) is given by 

S{iJ) ^ -^eilog2ei (14) 



(13) 



where e^'s are the eigenvalues of p{i,j). Using equation (9), (10), (11) (12) and (13), the first 
derivative of S{i,j) with respect to H is calculated near both the critical points (Fig. 5). The 
derivative diverges near Hc2 but approaches a finite value close to Hd- The n.n. concurrence can 
be written as [T4 l [34 1 135] 



a^^+i = 2 Max[0, \p23ii,i + 1)| - \/piiihi + l)/544(j, « + 1)] (15) 

Fig. 6 shows the derivative of Ci^i+i with respect to H versus H. The derivative, as in the case 
of one-site and two-site entanglement measures, diverges as H ^ Hc2 but has a finite value as 
H Hci . 

Lastly, we probe the existence of special features, if any, near the QCPs exhibited by the one-site 
RF [Ml ESI [301 EB Ea [33] defined in Eq. (5). The reduced fidelity susceptibility (RFS) is defined 
to be 

Xr[H) = hms^a (16) 

Figs. (7) and (8) show that the RF Tr{H, H + 5) drops sharply at the quantum critical points 
(insets) and the associated RFS, xb.{H), blows up as both the quantum critical points are ap- 
proached. This result is in contrast with what is observed in the case of entanglement measures, 
where a special feature develops only in the vicinity of the critical point Hc2 ■ The calculations of 
the RF and the RFS are possible because they involve only local measures. A calculation of the 
global fidelity would not have been possible lacking a knowledge of the true many body ground 
state. 
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FIG. 1: A two-chain ladder with rung and intra-chain nearest-neighbour exchange coupUngs of 
strengths J± and J|| respectively. The indices 1 and 2 label the two chains of the ladder. 
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FIG. 2: Plot of Al, the spin gap, versus from numerical diagonalization data of the ladder 
Hamiltonian (Eq. (1)), L being the number of rungs in the ladder, (inset) ground state fidelity 
versus magnetic field H for L = 3,4, ...16. 
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FIG. 3: Plot of versus H (using numerical diagonalization data); (inset) variation of ^ with 
H. 
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FIG. 5: Plot of '^^j^''^ versus H near H = Hc2', (inset) plot of versus H near Hci- 
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FIG. 6: Plot of ^"^'^"^^ versus H near H — Hc2', (inset) plot of '^'^'j^'*''" versus H near Hd. 
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FIG. 7: Plot of RFS Xr{H) versus H near H ^ Hci; (inset) plot of RF Fb.{H, H + 5) versus H 
near H = Hci . 
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FIG. 8: Plot of RFS XR^H) versus H near H = Hc2; (inset) plot of RF Fr{H, H + S) versus H 
near H — Hc2- 
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IV. DISCUSSIONS 



In this paper, we consider a spin- 5, two-chain AFM ladder in an external magnetic field. The 
ladder system is known to exhibit QPTs at two critical values, Hd and i/c2, of the magnetic field. 
The ladder has a rich quantum phase diagram with a gapless LL phase separating two gapped 
phases. Both the spin-disordered state {0 < H < Hd) and the fully polarized FM state {H > Hc2) 
constitute gapped phases. Using a bosonization technique, it has been shown [44] that the spin 
gaps vanish at the critical points and the spin-spin correlation functions become long-ranged. As 
suggested in [TB I flT j fTS]. a second order QPT is characterized by a discontinuous/divergent first 
derivative of an entanglement measure with respect to the tuning parameter. Our computations of 
the first derivatives of the entanglement measures S{i), S{i,j) and Ci^i+i show that these quantities 
diverge only a.s H ^ Hc2 but remain finite as the other critical point Hd is approached. As 
discussed in [16], the first derivatives of one or more elements of the reduced density matrix p{i,j) 
with respect to the tuning parameter are expected to diverge at the critical points. From Eqs. 
(11)-(13), one can verify that this is the case as H ^ Hd and Hc2 with the divergent contributions 
coming from pii{i,j) and p44{i,j). The theorem in [16] regarding the discontinuity /divergence of 
the first derivative of an entanglement measure at a critical point links the behaviour to that of the 
first derivative of one or more elements of p{i,j). This is so provided a set of conditions is satisfied. 
We find that one of these conditions (condition (b)) is violated in the case of the two-chain ladder 
as H — > Hd- This is easily illustrated for the entanglement measure Ci^i+i (Eq. (15)). The first 
derivative '^'"^'^^ involves terms containing the factor ■m{H) which leads to a cancellation 

of singularities as H ^ Hd (see Eq. (9)). This is contrary to condition (b) in [16] so that the 
theorem is no longer vaHd. The cancellation of singularities does not occur as H ^ Hc2 (see Eq. 
(10)) so that signals the occurence of a QPT. In the case of the single-site entanglement, 

S{i), similar arguments show that the cancellation of the singularity occurs as H —^ Hd- The 
square root singularities in magnetization (Eqs. (9) and (10)) are generic to other AFM systems 
with spin gap Hke the spin-1 chain in a magnetic field [lH |45l |46] . Thus, the behaviour reported in 
this paper may be a general feature of a class of gapped Id AFM systems. As shown in our paper, 
the measures RF and RFS yield appropriate signatures as both the critical points Hd and Hc2 are 
approached and thus appear to be better indicators of QPTs in the case of systems which violate 
one or more conditions of the theorem in [16] . 
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